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We propose an alternative mechanism based upon dark matter (DM) interpretation for anomalous 
peak signatures in cosmic ray measurements, assuming an extended dark sector with two DM species. 
This is contrasted with previous effort to explain various line-like cosmic-ray excesses in the context 
of DM models where the relevant DM candidate directly annihilates into Standard Model (SM) 
particles. The heavier DM is assumed to annihilate to an on-shell intermediate state. As the 
simplest choice, it decays directly into the lighter DM along with an unstable particle which in turn 
decays to a pair of SM states corresponding to the interesting cosmic anomaly. We show that a 
sharp continuum energy peak can be readily generated under the proposed DM scenario, depending 
on dark sector particle mass spectra. Remarkably, such a peak is robustly identified as half the mass 
of the unstable particle. Furthermore, other underlying mass parameters are analytically related to 
the shape of energy spectrum. We apply this idea to the two well-known line excesses in the cosmic 
photon spectrum: 130 GeV 7 -ray line and 3.5 keV X-ray line. Each observed peak spectrum is well- 
reproduced by theoretical expectation predicated upon our suggested mechanism, and moreover, 
our resulting best fits provide rather improved values. 


I. INTRODUCTION 

Dark matter (DM) is a necessary component in the 
time evolution of our Universe, and its existence has been 
supported by a tremendous amount of astrophysical and 
cosmological evidence for which the relevant observations 
are made mostly based upon the gravitational interaction 
of DM (see Ref. [T] for a general review). In fact, none 
of the Standard Model (SM) particles can explain vari¬ 
ous DM-related phenomena, so that the detection of any 
DM candidates can be not only exciting per se. but a 
strong sign of new physics framework beyond the Stan¬ 
dard Model (BSM). A great deal of experimental effort to 
detect DM signals has been made in three different direc¬ 
tions such as i) direct detection experiments by observ¬ 
ing recoil energy of target nuclei scattered off by DM, ii) 
indirect detection experiments by observing cosmic rays 
originating from DM annihilation or decay, and iii) col¬ 
lider searches (for example, at the Large Harden Collider 
at CERN) by actively producing DM candidates and ex¬ 
ploiting their collider signatures. These three avenues to 
DM detection are complementary to one another, and 
have set the bounds of the viable DM mass and the as¬ 
sociated cross section (see Ref. [5] and references therein 
for review).^ 

Of those experimental attempts, indirect detection ex¬ 


periments have received particular attention as many of 
them have reported anomalous observations potentially 
signalling the presence of DM candidates at the locus of 
cosmic ray sources. For instance, PAMELA [4], Fermi- 
LAT [S], and AMS-02 [0] found quite a marked rise of 
the positron fraction in the energy range from roughly 
10 to 200 GeV, and similarly ATIC [7], FERMI-LAT [8], 
and HESS [9] reported an excess in the positron-electron 
combined energy spectrum between 100 and 1000 GeV. 
Several photon channels showed intriguing excesses such 
as 3.5 keV line [101111], 511 keV line [12], GeV bump [13], 
and 130 GeV line [HI [15]. The positron excesses and the 
Galactic Center (GC) GeV 7-ray excess are featured by a 
continuum bump, while the other three X/7-ray excesses 
showed a sharp peak within a very narrow energy range. 

The latter class of excesses are particularly interest¬ 
ing because they can be readily connected to the DM 
interpretation. As typical DM candidates behave non- 
relativistically, the photon energy from a DM pair anni¬ 
hilation (or 2-body decay) is monochromatic, being the 
same as (half) the DM mass.^ In this context, many DM 
models to address those excesses have been introduced 
and studied in literature: for example. Ref. [H] for 511 
keV line. Refs. [T71 HH] for 130 GeV line, and Ref. [T5j 
for 3.5 keV line. In reality, the relevant signal spectrum 
does not appear as a J-function-like peak but is smeared 
to some extent because of imperfection in cosmic ray de¬ 
tectors. With the assumption of a Gaussian smearing. 


* email: immworry@ufl.edu 
^email: logl079@gmail.com 

^ Recently, Ref. proposed a general scenario, dubbed “Inflatable 
DM models”, within the context of which many well-motivated 
DM models having too large production of DM can be remedied, 
hence evade the bounds without tuning of underlying parame¬ 
ters. 


^ The 511 keV 7 -ray peak comes from the positronium decay. 
Thus, for the explanation of the 511 keV line excess, the required 
is a new source of positrons which can be DM annihilation or de¬ 
cay. 




the resultant 7 -ray energy spectrum (typically) becomes 
symmetric with respect to the nominal peak. 

This broadening effect has motivated the possibility 
of non-minimal DM scenarios for interpreting the nar¬ 
row width of peaks as the physical, not the one induced 
from instrumental uncertainties. The next-to-minimal 
DM models hypothesize the situation where DM parti¬ 
cles annihilate or decay into on-shell intermediate parti¬ 
cles that decay into photons nnun]. As such an on-shell 
intermediate particle comes along with a fixed boost, the 
photon energy spectrum is characterized by a rectangular 
shape [ 20 I [21] . If the mass gap between the DM and the 
intermediate particle is sufficiently small, hence so is the 
boost factor, then the photon energy spectrum becomes 
narrow enough, potentially being indistinguishable from 
the signal spectrum by the minimal scenario. 

Having similar philosophy and positing the DM inter¬ 
pretation, we here propose a new mechanism to develop 
a narrow continuum energy spectrum which would fake 
a sharp spike. The research program to explain the ex¬ 
cesses in cosmic ray energy spectra with the ‘‘‘'energy- 
peak" emerging under non-minimal DM frameworks has 
been initiated by Ref. [12], in which various observa¬ 
tions of the energy-peak made in the context of collider 
physics [23H28] have been applied to the GC 7 -ray GeV 
excess. As in [22], we begin the discussion with noting 
that multiple DM species could exist in the Universe, and 
the relevant DM models constructed upon such a DM 
framework can bring about not only nontrivial cosmo¬ 
logical implications, e.g., “assisted freeze-ouf [H], but 
interesting phenomenology, e.g., “boosted DM" [30H33- 
In this context, the proposed mechanism involves a non- 
minimal dark sector containing multiple DM particles. 

For the purpose of simplicity, we introduce two DM 
species one of which is assumed heavier than the other, 
denoting henceforth the former and the latter as Xh and 
Xl^ respectively. The heavier DM communicates to the 
SM sector not directly but through the lighter DM. In 
addition, the heavier DM pair-annihilates into a pair of 
on-shell intermediate states (denoted as A) each of which 
subsequently decays into the lighter DM together with a 
dark pion or an axion-like particle (ALP) (denoted as fU 
emitting a couple of photons in the final state.^ FIG. 
schematizes the DM scenario that we consider through¬ 
out this paper.^ We point out that although we employ 
the photon final state as a concrete example for elabo¬ 
rating our mechanism, it is straightforwardly extensible 
to other visible particle final states, e.g., e+e“. 

In this DM scenario, the on-shell intermediate particle 


® In general, A and xi can be either dark or SM sector parti- 
cles (and may be even unstable) unless they are stringently con¬ 
strained by other observations. However, for the sake of simplic¬ 
ity, we assume that they are dark sector particles. 

The main idea in this paper can be readily applied to the decay¬ 
ing DM scenario, but we keep the annihilating one as a concrete 
example. 



FIG. 1: The dark matter scenario under consideration. 

A comes with a fixed boost factor, leading to a rectan¬ 
gular energy spectrum, hence a rectangular boost dis¬ 
tribution for particle a. Due to variations in the boost 
factor of a the emitted photons develop a continuum en¬ 
ergy spectrum whose width is determined by the mass 
parameters involved in the entire process demonstrated 
in FIG. m We remark that the relevant scenario in¬ 
vokes from the (relatively) narrow energy spectrum to 
the broadly-distributed energy spectrum, depending on 
details of the DM models of interest. We shall briefly dis¬ 
cuss the aspect of the wide continuum bump signature in 
Sec. |III[ while primarily focusing on the “peak-faking” 
interpretation and the associated mass spectrum. 

Not surprisingly, the mass spectrum of ^ tua ^ 
ma+rn^i renders a relatively narrow photon energy spec¬ 
trum. As a characteristic feature, the resulting differ¬ 
ential energy spectrum is symmetric in the logarithmic 
scale, and remarkably, its center position is identified as 
half the mass of particle a, TOa/2 |23l [33]. Therefore, 
these structural properties enable us to not only distin¬ 
guish this DM model scenario from other standard DM 
interpretations, but probe/measure the mass parameters 
of some dark sector particles. 

To present our main idea, this paper is organized as 
follows. We first begin with the DM model under con¬ 
sideration in the next section. In Sec. |III[ we discuss the 
energy spectra of relevant visible particles arising from 
the DM model introduced in Sec. jll] mainly focusing on 
their functional structure. We then apply the main idea 
to a couple of photon peaks in Sec. |IV[ i) Fermi-LAT 130 
GeV excess and ii) 3.5 keV excess. Sec. is reserved for 
our conclusions. 


II. DARK MATTER MODEL 

We here discuss a class of DM models for which our 
main idea is applied, and offer a viable DM model that re¬ 
alizes the relevant scenario. As briefly mentioned earlier, 
we imagine that the dark sector is non-minimal, meaning 
that there exist more than one DM candidate. Although 
the arbitrary number of DM species could be introduced, 
we employ only two types of DM particles, the heavier 
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A. Derivation of the analytic expression 


(Xh) and the lighter (yz), for simplicity. 

The lighter DM is assumed to directly communicate 
to the SM sector, whereas the heavier DM is set out to 
have interactions with the SM sector via the lighter DM. 
In this sense, the relic abundance of Xh can be computed 
by the scheme of ^‘assisted freeze-ouf [H]. We further 
assume that Xh has a contact with xi directly, but 
via an on-shell intermediate state A, i.e., a pair of heav¬ 
ier DM particles annihilate into a pair of A’s each of 
which subsequently decays into xi together with a dark 
pion or an ALP as depicted in FIG. Finally, the dark 
pion or the ALP decays into a photon pair whose energy 
spectrum is the major interest of this paper. We remark 
that it is not a necessary condition for particle A to get 
pair-produced, i.e., A could be produced in association 
with another particle A' as the detailed dynamics of A! 
is irrelevant to the later argument and formalism. Sim¬ 
ilarly, A does not need to directly decay into the lighter 
DM, i.e., xi could be replaced by other heavy particles 
whose detailed dynamics does not affect the later argu¬ 
ment and formalism. In this context, the model set-up 
demonstrated in FIG. is a simplified version. 

A possible realization of the DM scenario at hand can 
be summarized as follows. We consider two fermionic 
DM species, Xh and x/, with an intermediate fermionic 
state V'A and a singlet pseudo-scalar a (e.g., a dark pion 
or an ALP as mentioned earlier). Then, the effective 
Lagrangian required for the DM scenario exhibited in 
FIG. [2 is simply described by the following operators: 

-Cdm a ^ XhXh'^Ai’A + A atpji'j^xi + Y (1) 


where F^’' denotes the dual field strength tensor as usual, 
and A and fa describe the associated suppression scales 
whose details can be revealed by appropriate UV com¬ 
pletion. The first term ensures an s-wave annihilation of 
the heavier DM, i.e., x^Xh -A ipA'f’A^ second induces 
the decay of if a into xi s-nd a, and the last corresponds 
to two photon decay of a as in FIG. respectively. The 
stability of Xh and xi can be easily achieved with sepa¬ 
rate symmetries, e.g., U(l)'G U(I)" [29] or Z 2 ® Z'f [^ . 
We again stress that the Lagrangian in Eq. 0 is a simple 
realization and there exist a host of other possibilities to 
accommodate the event topology in FIG. Exhausting 
all of them is, however, beyond the scope of this paper. 


III. ENERGY SPECTRUM 


Assuming that the heavier DM particles are non- 
relativistic, their pair annihilation into two A’s leads to a 
fixed boost of particle A (denoted as 7 ^) relating the two 
mass parameters by XA = /""^A with rui symbolizing 

the mass of particle species i. Since A obtains a non-zero 
boost factor, the energy of particle a is not monochro¬ 
matic, but given by a broad spectrum. The a energy. 
Eat measured in the laboratory frame is parameterized 
as 


E„ = E* 


lA ■ 


HE 

E: 


Y 72 -lcos6»; 


( 2 ) 


where 0* is the emission angle of a in the A rest frame 
with respect to the boost direction of A and A* is the 
fixed a energy measured in the rest frame of particle A: 


E* = 


m A — mr 


2mA 


(3) 


If A is either a scalar or produced in an unpolarized way, 
then cos 9* becomes a flat variable, resulting in a rect¬ 
angular distribution in Ea by a simple chain rule whose 
range is given by 


Ea € 


KIA - P*a \Jxa-^^ KIA + P*a \Jia-^ 


(4) 


Similarly to Eq. ([^, the observed photon energy for a 
fixed 7 q is expressed as 

E^ = E; {xa + ^Xl - 1 cos 0 ;) , (5) 

where 0 * denotes the intersecting angle between its emis¬ 
sion direction and the boost direction of particle a in the 
a rest frame and A* is the fixed photon energy measured 
in the a rest frame, that is, half the mass of particle a: 

E*, = ( 6 ) 

Unlike XA in the case of particle a, Xa is not single-valued 
but distributed. Denoting its distribution by g{xa)t from 
Eq. @ we find the unit-normalized expression (or equiv¬ 
alently, probability distribution function) for g{xa) as 

9iXa) = -Xa), (7) 

^P*aVlA-^ 

where 0(a:) is the usual Heaviside step function, and 7 ^ 
are defined by 


In this section, we derive the analytic expression for 
the gamma-ray energy spectrum and discuss its proper¬ 
ties which could be distinguished from other (standard) 
scenarios, given the scenario in FIG. 


X^^^XA±^Jxl-l- ( 8 ) 

ma ma V 

Here we used the fact that cos0* is a flat variable so 
that g{xa) develops a rectangular distribution as well. 
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We then find that for any fixed 7 a, the unit-normalized 
differential energy distribution is 


1 dr 

r dE^ 


fixed7a 




Q{E^ - E-)Q{E+ - E^)i9) 


where E^ can be obtained by setting cosd* to be ±1 in 
Eq. ([5|: 

E^ = E;{ja ± - !)• ( 10 ) 


Denoted by f{E^), the expression for the unit-normalized 
total energy spectrum can be obtained by summing 
Eq. over all relevant 7 a’s, that is, 


fiEy) = f 


J aila) 

d'la - , 

7 S"” - 1 


( 11 ) 




,W 2 _ V(7““)^ - 1 + 7a 

'avlA 


-log y/(7““)^ - 1 + }. 


( 12 ) 


where 7 ™™ and 7 “®'’' are defined as 


7 ““ = max 


7a 


1 fE.E; 


2 \E* E^ 


rr"^7a+- (13) 


Since 5 ( 70 ) is upper-bounded, 7 ^ determines the span¬ 
ning range of f{Ej) as follows: 


Ey 

— e 
Et 


it - \/ilt)"^ - 1 , it + \/(7^)^ - 1 


(14) 


We finally remark that in the actual data analysis with 
concrete examples, all prefactors in Eq. are even¬ 
tually absorbed into the overall normalization parameter 
N, and as a consequence the shape of f{E^) is completely 
determined by 7 +, 7 “, and E*, i.e., there are four inde¬ 
pendent fit parameters. 


B. Functional properties and discussions 


To discuss functional properties of the photon energy 
spectrum, we first revisit the expression of Ej for a fixed 
7 a shown in Eq. (j^. Since cos 9* spans —1 to -|-1, the 
range of E-y is trivially given by 


E^ 

El 


la - \/ll-^. la + \/ll - 1 


(15) 


as also expressed in Eq. (10). One remarkable feature 


from the above range is the fact that the lower (upper) 
end in the right-hand side is smaller (greater) than 1 , 
implying that E* is the only commonly-included energy 
value for any 7 a. Moreover, we observe that E* is the ge¬ 
ometric mean of minimum and maximum energy values, 
i.e., {E*Y = This again implies that the E* 


value is located at the center of the E^ distribution for a 
given 7 a in the logarithmic scale. 

As mentioned before, the shape of E^ distribution for 
a fixed 7 a is rectangular due to the flatness of cos 9* in 
Eq. (§. In order to obtain the overall energy distribu¬ 
tion, one should “stack up” all of such rectangles con¬ 
tributing to the energy distribution. This statement is 
already formulated in Eq. 0 as a Lebesque-type inte¬ 
gral representation. Hence, one would expect that the 
resultant photon energy distribution contains a unique 
peak at E^ = E* in conjunction with the observation 
made in the previous paragraph. Indeed, there arises a 
subtlety here: the validity of this expectation depends 
whether or not the smallest boost factor of particle a 
approaches 1 . 

The condition for 7 a = 1 (i.e., solving Eq. ([^ with 
7 “ being 1 ) is A* = rria'iA^ representing a hypersur¬ 
face formed by rriA, and for any fixed rua, i-S., 
only delicately chosen mass spectra can attain this con¬ 
dition. Once 7 a = 1 is available, clearly, E* appears as a 
unique and cusp-structured peak [53] . The left panel 

in FIG. [^demonstrates an example spectrum in this cate¬ 
gory in the logarithmic scale. The chosen mass spectrum 
is niA, rua) = (237.5, 200, 50, 100) GeV, 

and events were generated by pure phase space. The sim¬ 
ulated events are binned to the blue histogram, and the 
associated theory prediction is shown by the red line. We 
clearly see that the theory expectation can reproduce the 
data well enough, and the spectrum is symmetric with 
respect to E* = ma/2 = 50 GeV (indicated by a black 
dashed line) in this scale. We also remark that both sides 
of the distribution loo k li ke straight lines, which can be 

in 

+ E-,) in 

he whole spectrum appears as 


easily seen from Eq. (12) with 7 “''' = 
logarithmic E^^ so thaf 
an isosceles triangle. 

On the other hand, g{"fa) starts from the value away 
from 7 a = 1, 710 peak is developed in the middle of the 
energy distribution. Instead, a plateau region emerges 
because even the narrowest rectangle corresponding to 
the smallest 7 a has a finite-sized width. Nevertheless, it 
is straightforward that the center of the relevant photon 
energy spectrum can be identified as E* in the logarith¬ 
mic base. These expectations are manifestly shown in 
the right panel of FIG. [^ The relevant mass spectrum is 
the same as the previous case with mA replaced by 170 
GeV. Also, the existence of a plateau region makes the 
whole energy spectrum appears as an isosceles trapezoid 
in the logarithmic base. This plateau structure is a dis¬ 
tinguished feature from other energy spectra. However, 
its presence may be invisible in actual indirect detec¬ 
tion experiments due to the issue of energy resolution. 
Clearly, if the plateau is smaller than the relevant re¬ 
solving power, its existence is rarely identifiable so that 
the relevant spectrum easily fakes a unimodal distribu¬ 
tion like the previous case. Even for the energy spectrum 
with a sufficiently sizable plateau, its identification may 
be unavailable with small statistics, i.e., more data accu¬ 
mulation may be required. 
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Case with gi;'a)tO 



Case with a plateau 



[GeV] 


FIG. 2: Left panel: the gamma-ray energy spectrum with a peak. The chosen mass spectrum is rriA, , rna) = 

(237.5, 200, 50, 100) GeV. The simulated data and corresponding theory expectation are represented by the blue histogram 
and red line, respectively. Right panel: the gamma-ray energy spectrum with a plateau. The chosen mass spectrum is the same 
as in the left panel with ruA replaced by 170 GeV. 


As the proposed mechanism, which is denoted as Sce¬ 
nario iii), is used for explaining narrow peaks, it is in¬ 
teresting to compare and contrast it with the other two 
conventional scenarios enumerated below. 


• Scenario i): Photons directly come from DM anni¬ 
hilation/decay. 

• Scenario ii): Photons are emitted as a decay prod¬ 
uct of the on-shell intermediate particle into which 
DM annihilates/decays. 


In Scenario i), the width of the peak is typically caused by 
the intrinsic energy resolution of detectors, and thus the 
final energy spectrum is symmetric about the nominal 
peak. Even in the logarithmic scale, the spectrum is de¬ 
scribed by a smooth curve, while it is no more symmetric- 
looking about the peak. In Scenario ii), the width of the 
energy spectrum is physical. However, identifying the 
peak position is ambiguous due to the box-like spectral 
behavior in both the linear and the logarithmic scales. 
Therefore, defining the symmetry property of the shape 
is not available. The comparisons thus far are summa¬ 
rized in Table [ij One can easily see that their respec¬ 
tive morphological features differ from one another, and 
therefore, one is able to pin down the underlying DM 
scenario with a reasonable amount of signal statistics. 


Before closing the current section, we briefly discuss 
the case of wide continuum bump spectra. Certainly, 
such types of spectra arise within a broad realm of rel¬ 
evant parameter space so that the theory prediction in 
Eq. (12) can be employed to explain continuum cosmic 
ray excesses. Due to the unique morphological features 
discussed so far, the relevant signal spectrum could be 
easily distinguished from other continuum bumps, taken 
as strong evidence of a non-trivial dark sector. 



Scenario i) 

Scenario ii) 

Scenario iii) 

Peak existence 

Always 

Absent 

Sometimes 

Plateau existence 

Absent 

Always 

Sometimes 

Width 

Instrumental 

Physical 

Physical 

Symmetry in E 

Symmetric 

Not available 

Asymmetric 

Symmetry in log A 

Asymmetric 

Not available 

Symmetric 

Shape in E 

Gurved 

Rectangular 

Curved 

Shape in log E 

Gurved 

Rectangular 

Oblique 


TABLE I: Comparisons of structural properties in the en¬ 
ergy spectrum among three DM scenarios defined in the text. 
Symmetry properties are defined with respect to the relevant 
peak if available. 


IV. APPLICATIONS 

Armed with the argument in the previous section, we 
apply the basic idea to a couple of existent cosmic ray 
peaks: i) 130 GeV line [HKTS] and ii) 3.5 keV line mm 
Although claiming that these examples could be under¬ 
stood by the proposed mechanism, we admit that the 
underlying DM models for them might not fall into the 
scenario depicted in FIG. Nevertheless, we empha¬ 
size that the applicability of the relevant technique is 
restricted to neither these two examples nor gamma-ray 
spectra. In other words, nothing precludes us from ap¬ 
plying the DM interpretation at hand for any of future 
excessive cosmic ray signals with a narrow peak. 


A. 130 GeV line 

Our first example is the famous 130 GeV line whose 
original data was collected by the FERMI-LAT collab¬ 
oration. We basically conduct the fit to the spectrum 
of the observed 130 GeV gamma-ray excess with the ex- 
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pected shape in Eq. with the overall normalization From Eq. ® , we obtain the expressions for E* and p* in 
parameter Ns being added: terms of 7 ^ 


fsiE^) = iVsjlog V ( 7 ““)^ - 1 + 7r 

-log [v/(7 r")"-l + 7“”]}- (16) 


The relevant data points and errors are taken from Reg4 
with the ULTRACLEAN event class in Ref. [U]. Since 
the measured bin counts contain the contributions from 
backgrounds, the relevant fit is performed simultaneously 
with the background template. We assume that the back¬ 
grounds can be parameterized by a simple (gradually- 
falling) power law such as 




(17) 


where Nb is the normalization parameter for back¬ 
grounds and p encodes the background shape. Here E* 
is the same E* contained in 7 ™“ of Eq. (16). There¬ 
fore, the entire data set is fitted with the combination of 
fs{E^) in Eq. (fl^ and fsiE^) in Eq. (1T7|, i.e.. 


ftoUEj) = fs{E^) + fsiE^). 


(18) 


Our fit result is shown in the upper-left panel of FIG. 
The data points and associated error bars are represented 
by black dots and black lines, correspondingly, while the 
red solid curve shows the best-fitted model. We clearly 
see that our fit is in a rather good agreement with the 
gamma-ray energy spectrum, i.e., the template given in 
Eq. (18) reproduces the data sufficiently well. More 
quantitatively speaking, we find that the value is 19.7 
for 19 degrees of freedom (i.e., 25 data points subtracted 
by 6 fitting parameters such as Ns, Nb, 7 ^, 'Ja , E*, and 
p) between 80 and 210 GeV. This number is quite com¬ 
parable to that from conventional scenarios, suggesting 
that our DM scenario be considered (equally) plausible 
in explaining the observed data. 

The fit also tells us useful information. First of all, we 
extract the mass of particle a from the measurement of 
E*{= mal2), that is. 


m-- = 258.1Ii 4';2 GeV, 


(19) 


where the errors here are evaluated by Icr statistical un¬ 
certainty. The other mass parameters can be also esti¬ 
mated by the measured 7 ^ and 7 “, for which the best- 
fitted numbers are reported as 


,1-1-0.0418 


(7a Qg4g, 

( 7 -)®='* = 1 . 0000 +°•°°^^ 


( 20 ) 

( 21 ) 


respectively. Again, the errors are reported by Icr statis¬ 
tical uncertainty. Note that the lower error for ( 7 “is 
not provided because there is no sensitivity to the v alue s 
of 7 - < 1 according to the definition of 7 “'" in Eq. (131. 


EllA 



»^a(7a +7a ) 

2 

mailt -la) 

2 


( 22 ) 

(23) 


and the difference between the former squared and the 
latter squared leads the following mass relation: 


o -L “T 


2mAma 


= ltla 1 


(24) 


with which we perform a scan of allowed mass space. 

As 7 - shows the smallest error, we simply fix it to 
be unity for phenomenological purpose. Instead of do¬ 
ing three-dimensional scanning, we choose three different 
values, 10, 100, and 1000 GeV to find the allowed 
regions in terms of and niA- The lower-left panel of 
FIG. [^demonstrates the allowed parameter space in the 
plane of {niA — — nia ) versus — ttia — — iria) ■ 

The red, green, and blue regions are for = 10, 100, 
and 1000 GeV, respectively. Solid curves are the contours 
evaluated with the best-fitted nia and 7 +. On the other 
hand, dashed (dot-dashed) curves are the ones evaluated 
with nia + Snia and 7 + -I- ^ 7 + (nia - Snia and 7 + - 1 ^ 7 +), 
so that one can get some intuition on the mass spectrum 
allowed by Icr variations of the relevant parameters. We 
observe that the viable mass spectra are more or less 
compact. This is not surprising because the narrow peak 
enforces a degenerate mass spectrum not to obtain too 
large boost in any of the steps in FIG. Dark sec¬ 
tor scenarios featuring such a mass spectrum could be 
achieved by a symmetry. Building realistic DM models 
is, however, beyond the scope of this paper, so we do not 
further pursue this direction here. 


B. 3.5 keV line 


Our next example is the well-known 3.5 keV line. As 
in the previous case, the fit is performed to the spectrum 
of the observed 3.5 keV gamma-ray excess with the sig¬ 
nal template given in Eq. (16). The relevant data points 


and associated errors are taken from the processed data 
for the MOS spectrum of the central region of the An¬ 
dromeda galaxy (M31) found in Ref. [TT]. The numbers 
to be used contain only the DM component, i.e., the back¬ 
ground component is already subtracted. Therefore, we 
execute the fit procedure only with the signal template, 
not introducing any background template unlike the pre¬ 
vious case. 

The fit result is exhibited in the upper-right panel of 
FIG.[^ Again, the data points and the error bars are rep¬ 
resented by black dots and black lines, respectively, while 
the red solid line describes the best-fitted model. Our fit 
is in a very good agreement with the measured energy 
spectrum. The reported value is 4.61 for 13 degrees 
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FIG. 3: Upper-left panel: the 130 GeV 7 -ray energy spectrum of the DM and background components taken from Ref. m- 
The fit is performed with 25 data points (black dots), and the best-fitted is represented by the red curve. Upper-right panel: 
the 3.5 keV X-ray energy spectrum of the DM components taken from Ref. m- The fit is performed with 17 data points. 
Lower-left panel: the allowed mass space for the 130 GeV line along with la variations of ma and 7 ^ for three different masses 
of particle Xi- Lower-right panel: the allowed mass space for the 3.5 keV line. 


of freedom (i.e., 17 data points subtracted by 4 fitting 
parameters such as Ns, la ^ E*) between 3 and 

4 keV. This number is significantly improved, compared 
with that from standard interpretations, and therefore, 
our DM framework can be considered as a plausible sce¬ 
nario accommodating the observed spectrum. 

Speaking of various best-fit parameters, we first find 
that the extracted mass parameter for particle a is 

= T.OQtHe keV, (25) 

where errors are estimated by Icr statistical uncertainty. 
We find the best-fit values for 7 + and 7 “ are 

( 7 ,+ )'^’^* = (26) 
( 7 -)®='* = 1.0000+°•°°^^ (27) 

respectively, together with Icr statistical uncertainty. In 


order to obtain the allowed mass space within Icr vari¬ 
ations of the relevant parameters, we follow the same 
procedure in the previous case for three different 
masses, 1, 10, and 100 keV. The lower-right panel of 
FIG. [^demonstrates the allowed region again in the plane 
of {mA — — ma) versus (m^,^ — mA — — ma)- The 

red, green, and blue regions correspond to m^^ = 1 , 10 , 
and 100 keV, respectively. 


V. CONCLUSIONS 

Dark matter indirect detection offers an excellent op¬ 
portunity to confirm the existence of DM. Several exper¬ 
imental collaborations have already reported anomalous 
phenomena in the relevant cosmic ray energy spectrum. 
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Particular attention has been paid to sharply-peaked sig¬ 
nals due to the easiness of DM interpretation. Typical 
models assume that (non-relativistic) DM particles di¬ 
rectly annihilate or decay into visible ones. In this min¬ 
imal DM scenario, the narrow width of the peak is typ¬ 
ically understood as the one stemming from the imper¬ 
fection of cosmic ray detectors. 

We have rather taken the viewpoint that this width 
can be physically induced, and proposed a new mecha¬ 
nism to realize it with the assumption of a non-minimal 
dark sector. Two DM species were introduced, and the 
heavier one is assumed to annihilate to the on-shell inter¬ 
mediate state which subsequently decays into the lighter 
one and an unstable particle. This unstable particle fur¬ 
ther decays into a pair of visible particles which can be a 
source of anomalous peaks in the cosmic ray energy spec¬ 
tra. We showed that the signal spectrum can be narrow 
enough to fake a sharp spike with a suitable choice of the 
associated mass parameters. 

The shape of the full signal spectrum was derived, and 
several interesting functional properties were discussed. 
We pointed out that the peak position, one of the fit pa¬ 
rameters, is immediately identified as half the mass of 
the above-mentioned unstable particle. We also showed 
that other mass parameters can be estimated by other fit 
parameters. We then enumerated various morphological 
features to be utilized for distinguishing several DM sce¬ 


narios in which the sharp peak signature is available. The 
viability of the relevant strategy was assessed with two 
real observational data sets, 130 GeV line and 3.5 keV 
line. We found that both of them can be well-described 
by the theoretical expectation in our DM scenario, and 
that each allowed parameter space is fairly large. 

Finally, we emphasize that our DM scenario and the 
associated data analysis are not restricted to photon en¬ 
ergy peaks, i.e., any cosmic ray energy peaks can take 
advantage of them. We strongly encourage people to pur¬ 
sue the direction exploited in this paper as well whenever 
cosmic ray peaks are observed. 
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